TPJU-15/02, IFUP-TH/2002-36, MPI-PhT/2002-45 



Exact Witten Index in D=2 supersymmetric 
Yang- Mills quantum mechanics 

§ ; M. Campostrini, 

^ I INFN, Sezione di Pisa, and 

D , Dipartimento di Fisica "Enrico Fermi" dell'Universita di Pisa, 



> 
o 



Via Buonarroti 2, 1-56125 Pisa, Italy, 
J. Wosiek, 

M. Smoluchowski Institute of Physics, Jagellonian University, 



^ ■ Reymonta 4, 30-059 Krakow, Poland 

O ■ September 17, 2002 

(N ■ 
O 

Abstract 

Q I A new, recursive method of calculating matrix elements of polyno- 

' mial hamiltonians is proposed. It is particularly suitable for the recent 

algebraic studies of the supersymmetric Yang-Mills quantum mechan- 
^ I ics in any dimensions. For the D = 2 system with the SU(2) gauge 

^ [ group, considered here, the technique gives exact, closed expressions 

for arbitrary matrix elements of the hamiltonian and of the supersym- 
metric charge, in the occupation number representation. Subsequent 
numerical diagonalization provides spectrum and restricted Witten in- 
dex of the system with very high precision (taking into account up to 
10^ quanta). 

Independently, the exact value of the restricted Witten index is 
derived analytically for the first time. 

PACS: ll.lO.Kk, 04.60.Kz 
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1 Introduction 



Supersymmetric quantum mechanics provides a simple laboratory to study 
many properties of supersymmetric systems [0, . Recently supersymmet- 
ric Yang-Mills quantum mechanics (SYMQM) in d space dimensions attracts 
a lot of attention because of its possible relation with M-theory Q for = 9 
and large number of colours. Even though the three loop calculations ^ 
question the exact equivalence between the two, it still remains a valuable 
model of the latter sharing many of its features (e.g., continuum spectrum 
of scattering states with the threshold bound state - a supergraviton) . At 
the same time, studies of the lower-dimensional systems with various gauge 
groups provide the global understanding of the whole family of these models 
with many interesting limiting cases. A good example is the D = A SYMQM 
with the SU(2) gauge group whose spectrum, in the zero fermion sector, is 
identical with that of the well known 0- volume glueballs 0, ||, 0] • More- 
over, supersymmetry guarantees that, in addition to the continuum of the 
scattering states, there must also exist localized gluino-glueball bound states 
with the same masses. Indeed such solutions were found in Ref. |T0[ where 
a new algebraic approach to study these systems was proposed. Compact 
version of the D = A SYMQM was also studied recently by van Baal who 
pushed rather far analytical understanding of this system [11], O]. 



The technique of IT^] is the adaptation of the hamiltonian methods ||T3 
to supersymmetric systems with local gauge invariance. It consists of two 
steps: first, a finite (e.g., cut off) basis of gauge-invariant eigenstates of the 
bosonic and fermionic occupation number operators B = atal and F = fbfl 
is generated (we omit all but colour indices for a moment). Second, the ma- 
trix representation of the hamiltonian (and any other relevant observable) is 
derived and numerically diagonalized. All this is automated by implementing 
standard rules of quantum mechanics in an algebraic language like Mathe- 



matica. A faster, compiler based, version is now available |T^. As the cutoff 
A'cut we choose the number of bosons B, i.e., the cut off basis consists of all 
states with B < N^nt and all allowed fermions. The cutoff is gauge and rota- 
tionally invariant and consequently the spectrum reveals the full SO(D— 1) 
symmetry for each Ncnt- The technique has been applied to Wess-Zumino 
quantum mechanics, D = 2 and D = A SYMQM and to D = 5-10 YMQM, 
all based on the SU(2) gauge group [|10], |T5|. In all cases studied until now the 
spectrum of lower states converges with iVcut before the number of states be- 
comes unmanageable. Calculations become more and more time consuming 
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with increasing D and but the D = 10 case is within the reach of today's 
computers for the first few values of A^. Similar methods have been indepen- 
dently developed in Refs p!6| , p!7[] to study lower dimensional supersymmetric 
field theories. 

In the first part of this letter we present a new method to calculate matrix 
representations and apply it to the D = 2 SYMQM in Section 3. We calcu- 
late recursively the matrix elements of the hamiltonian thereby eliminating 
lengthy and space consuming process of generation and storing of the basis. 
For the D = 2 system recursions can be solved resulting in closed expressions 
for any matrix element of the hamiltonian H. We then diagonalize H numer- 
ically and calculate the restricted Witten index for this system to virtually 
arbitrary precision. The method can be generalized to higher D, allowing to 
reach values of the cutoff larger than those in Ref. . 

The second part contains an exact calculation of the restricted index, 
exploiting the analytic properties of the densities of bosonic and fermionic 
states, suitably regularized in the infrared. To our knowledge this is the first 
analytic calculation of the restricted index for this system. 



Supersymmetric Yang-Mills quantum me- 
chanics in two dimensions 



This model is exactly soluble g] also for higher gauge groups [|T8[ . It is how- 
ever still interesting as it shares some of the complexity of higher dimensional 
models. For example it has a continuous spectrum which is a characteristic 
feature of many other supersymmetric models and poses some challenge for 
the hamiltonian methods. Moreover it has a nontrivial Witten index which 
was defined only recently [IlOi . 



The system, reduced from D = 2 to one (time) dimension |]T9|], is de- 



scribed by the three real bosonic variables Xa{t) and three complex, fermionic 
degrees of freedom ipait), both in the adjoint representation of SU(2) , a = 
1,2,3. 

The hamiltonian reads 

H = ^PaPa + igeabc'lplxbi^c, (l) 

where the quantum operators x,p, ip, 4''^ satisfy the canonical (anti) commutation 
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rules 

[Xa,Pb\ = iSab, {i'a^i^l} = Sab, (2) 

and can be written in terms of the creation and annihilation operators 

Xa = ^{aa + ai), pa = -^{aa-a}^), (3) 

= /a, = fl (4) 

The system has a gauge invariance with generators 

Ga = eabc{xbPc - ii^lipc)- (5) 

Therefore the physical Hilbert space consists only of the gauge-invariant 
states. This constraint is easily accommodated by constructing all possi- 
ble combinations of creation operators (creators) invariant under SU(2), and 
using them to generate a complete gauge-invariant basis of states. There are 
four lower order creators: 

{aa)^aial {af)^a\fl (af f) ^ eabcalflfl [f f f) ^ ^abcflfUl (6) 

Fermionic creators satisfy (//) = {afY = {affY = {f ffY = 0, therefore 
the whole basis can be conveniently organized into the four towers of states, 
each tower beginning with one of the following states 

|0^) = |0); |l^) = (a/)|0), |2^) = (a//)|0), |3^) = (///)|0), (7) 

where we have labeled the states by the gauge-invariant fermionic number 
F = fafa- To obtain the whole basis it is now sufficient to repeatedly act on 
the four vectors (|^) with the bosonic creator (aa). Acting with other creators 
either gives zero, due to the Pauli principle, or produces a state from another 
tower, already obtained by application of (aa). The basis with cutoff Ncnt is 
then obtained by applying (aa) up to Ncnt times to each of the four "base" 
states of Eq. (^. Obviously our cutoff is gauge-invariant, since it is defined 
in terms of the gauge-invariant creators. 

The hamiltonian (0) reduces in the physical basis to that of a free bosonic 
particle 

H = -PaPa + gXaGa, (8) 

therefore it preserves the fermionic number and can be diagonalized indepen- 
dently in each sector spanned by the four towers in Eq. (^. 
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The spectrum is doubly degenerate even at finite Ncnt, because of tlie 
particle-liole symmetry wliicli is preserved by our cutoff. Particle-liole sym- 
metry relates empty and filled fermionic states (|0f) ^ |3f)) and their 1- 
particle 1-hole counterparts ^ |2f))- On the other hand, the super- 

symmetry generator 

Q = T.^aPa, (9) 

a 

connects sectors which differ by 1 in the fermionic number, e.g., it connects 
Of — If sectors and 2f — Sf sectors, but it does not connect If and 2f sectors. 
This can be easily seen in terms of the "-B-parity", Pb = (—1)'^: Q not only 
changes F, but also Pb', however, the F = 1 and F = 2 sectors have the 
same Pb, hence matrix elements of Q between these sectors vanish. 
Witten index vs restricted index. Because the particle-hole symmetry inter- 
changes odd and even fermionic numbers, the Witten index 

Iw{T) = exp i-TE,), (10) 

vanishes identically for this model|5 Nevertheless one can obtain a nontrivial 
and interesting information by defining the index restricted to a one pair 
(e.g.. Of — If) of sectors 

J(T)(°'i) = S,,^^=o,i(-l)^' exp i-TE,). (11) 

Since Q does not connect the If sector with the 2f sector, supersymmetry 
balances independently, and identically, fermionic and bosonic states within 
the Of — If and 2f — Sf pairs, with the usual exception of the vacuum. There- 
fore the restricted index is a good and nontrivial measure of the amount of the 
violation of SUSY, even when the total Witten index vanishes. Obviously, 
j(2,3) _ _j(o,i) jj^ principle, one could also consider and however, 
similarly to the global index, they vanish due to the particle-hole symmetry. 
Studying the restricted index is particularly interesting in this model since, 
due to the continuum spectrum, it does not have to be an integer. 

3 Exact matrix elements and numerical diag- 
onalization 

In order to simplify numerical computations, it is useful to avoid dealing 
explicitly with the Fock space vectors. This can be achieved by writing every 

^This is also true for any finite cutoff since A^cut preserves tlie particle-liole symmetry. 
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quantity of interest as a vacuum expectation value (vev) of a gauge-invariant 
operator, and deriving recursive relations between such vev's. 

We shall deal here explicitly with the F = and F = 1 sectors, the F = 2 
and F = 3 sectors can be obtained by the particle-hole symmetry. To begin 
with, let us introduce the gauge-invariant operators 

^ = Y1 ^bttb, = Y1 ^Wb = ('^«)' ^ = = X! ^bol, (12) 

b b b 

which satisfy the commutation relations 

[A,A^=4B-6, [A,B] = 2A. (13) 
States in the F = sector have the form 

|2n,0) = ^(At)"|0), (14) 



and those with F = 1 can be written as 

|2n+l,l) = ^Ft(At)"|0), F = Y.abfb, Ft = ^ aj/t = (a/). 

V Cfi b b 

(15) 

Observe that in the F = sector fafl = 6ab and therefore FF'^ = B. We are 
interested in the scalar products 

(0|A"'(At)-|0) = 5„,„c„, (16) 
{0\A^'FF\A^r\0) = {0\A^'B{A^r\0) = 5n'nC (17) 

and in the matrix elements of H and Q. Since 

H = -\{A + A^ -2B + 3), (18) 

the matrix elements of H can be written in terms of the above vev's and of 
(0|A"'S2(At)"|0). Given Eqs. (D and (|), the matrix elements of Q are also 
expressible by the vev's defined in Eqs. (|1^) and ([T7|). 

A general technique to compute the desired vev's is to consider a generic 
matrix element 

(0|A"ieA"2H(v4t)"|0) 

= (o|A"^-^eA"2+^H(At)"|o) + (o|A"i-^[A,e]A"2S(v4t)"|o), (i9) 
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where O = B, and E = 1, B. Eq. (|T|) allows to "shift to the left" 6 until 
it is immediately to the right of (0|. Then we use (01^4"'' = or {0\B = (0|3; 
the commutator terms are of lower order and therefore the iteration of Eq. 
([T9| ) closes, giving the desired recursion. 

This method can be applied to SYMQM in arbitrary dimension. For 
D = 2, recursions can be solved, providing explicit expressions for the matrix 
elements of H and Q between orthonormalized states. To this end we consider 
first the F = sector and define 

(0|AMt = ln{0\A''-\ (20) 

In can be computed recursively: /i = 6 and 

/„(0|A"-i = {OlA^'-^AA^ = {0\A''-\A^ A + 4B - 6) 

= miln-i - 6)A"-i + 4(K-\ B] + BA^'-')] = {In-i + 8n- 2)(0| 

(21) 

therefore = 2n + 4n^. Moreover, c„ = /„ c„_i and Cq = 1. We next compute 

{2n,0\B\2n, 0) = — (0| A"E(A"^)"|0) = — (0|(2n + 5)A"(At)"|0) = 2n + 3 

(22) 

and 

{2n-2,0\A\2n,0) = {2n,0\A^\2n-2,0) = J = V2n + An^. (23) 

V Cn-1 

Therefore, the nonzero matrix elements of H are 

(2n, 0\H\2n-2, 0) = (2n-2, 0\H\2n, 0) = -iV2n + 4n2, (24) 
{2n,0\H\2n,0) =n + l. (25) 

The F = 1 sector is dealt with in the same way, and the result is 

(2n+l, l\H\2n-l, 1) = (2n-l, l\H\2n+l, 1) = -IVQn + An^ (26) 
{2n+l,l\H\2n+l,l) = n + l (27) 

For the supersymmetry generator we obtain 

(2m, 0|Q|2n+l, 1) = -i-y/m+ | + (28) 
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In the cut Hilbert space, with states up to the |2Ai'cut+2, 0) and \2Ncut+'^, 1), 
Q is represented by the (A^cut+2) x (A^cut+1) matrix Qij = {2i,0\Q\2j+l,l) 
and by its adjoint. Define 



(2^,0|i/cut|2j, 0) = '^J2Q^kQlj, (2z+l,l|ifeut|2j+l,l) = lJ2QlkQkj- 

k 

(29) 



k k 

All matrix elements of H and iJcut are equal, apart from 

(2iV,,t+2, 0|/7cut|2iVe,t+2, 0) = \{N,^, + 1), (30) 
cf . Eq. (^) . Hence we confirm that we can define a cut off hamiltonian if^ut 



whose spectrum has exact supersymmetry at any finite A^cut llOj 0, p^ - 

The hamiltonian matrix has a tridiagonal structure and it can be diago- 
nalized numerically in a very efficient way. We use the 0{N'^) algorithm im- 
plemented in the lapack library, which computes all eigenvalues for N^^t = 10^ 
in a few minutes on a PC. 

A plot of the restricted index as a function of the Euclidean time is 
presented in Fig. |I] confirming results of |TD| to a much higher precision, and 
strongly suggesting the exact, time independent, value I{T) = 1/2. 

Note the intriguing multiple crossing at T = 1 where the index seems to 
attain its asymptotic value for all A^cut • Upon closer inspection it turns out 
that the curves do not cross at the same point. However the intersection 
points approaches Iw = 1/2, T = 1 as 6Iw = 0{1/Ncut), = 0(l/\/Acut)- 
This suggests existence of a "duality" transformation which relates, for finite 
A^cut, SUSY violations above and below the "critical" energy E ^ 1. 



4 Exact calculations of the restricted index 

Gauge-invariant eigenstates of the free hamiltonian, Eq. (^ are labeled by 
the absolute value of the momentum p. In this representation the index 
readsQ 

J(T)= / {pB{p)-pF{p))e-^^'^'dp, (31) 
Jo 

Since the spectrum is continuous, definition of the bosonic and fermionic 
densities, pb and pp, requires an infrared regularization. Therefore we first 

^We consider the index restricted to the F ~ and F = 1 sectors but will omit for the 
simplicity the (0, 1) superscript. 
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Figure 1: The restricted Witten index computed for values of Ncut ranging 
from 125 to 128000. 



solve the free problem in a spherical well of radius R and then calculate the 
index as a limit^. 

J(T) = hm /«(T), (32) 
with the IR regularized index given by the discrete sum 

/^(T) = Y: (e-^(^")'/2 - e-^(P")V2) , (33) 

n 

where p^{p^) are the discrete momenta of the bosonic and fermionic states 
respectively. 

The spectrum of a free particle in a three dimensional spherical well of 
radius R is determined by the boundary condition ji{pR) = 0, for the l-th 
spherical wave. However the local gauge symmetry limits allowed values of 
/ to / = 0, 1. Only these two angular momenta can be combined with the 

■^We thank Ken Konishi for the discussion on this point. 
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fermionic colour spin, s = 1, into a scalar, J = 0, gauge-invariant state in 
the F = and F = 1 sectors.^ It follows that the sum in Eq. (|33|) is over the 
positive zeroes z![^ of the first two Bessel functions. 



B 



(1) 



A simple way to compute the R ^ oo limit of Eq. 
by its large-n asymptotic expansion 



(34) 

is to replace z^^ 



7rn, 



zi^) =(3n- l//3n + 0{n-'), I3n = 7i{n + i); 



(35) 



we note that that, in the R —>■ oo limit, higher-order terms do not contribute 
and the sum can be replaced by an integral: 



j dn |exp 


2 2 

—u n 


— exp 


" 2^2' 


exp 


'-u\n + \f 










7r2 _ 







TT 



lim si — exp 

u^o 2u 1 



2u^ 



71^ 



erfc 



1 

2' 



(36) 



where 



u 



2^' 



2 f°° 

erfc(z) = / exp[— t ] dt, erfc(O) = 1. 
\/7r Jz 



(37) 



A more illuminating strategy to compute Jw starts with the well known 
theorem on meromorphic functions to rewrite Eq. ( p3D as a contour integral 



iRiT) 



1 

27ri Jr 



-T{z/R)y2 



h{z)dz, 



(38) 



where F denotes any contour enclosing counterclockwise positive real axis 
and contained within the angular region —tt/4 < arg(2) < 7r/4. The weight 
h is given by 

h{z) =j',{z)/Uz)-j[{z)/Mz) + II z. (39) 

The first two terms follow directly from the above theorem, the last one 
subtracts the pole at ^ = 0. Since ji{z) has a second-order zero at 2; = 
the theorem would not apply. However we can subtract by hand the singular 



*The "angular momentum" considered here is the colour angular momentum which 
generates rotations in the three dimensional colour space. 
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terms of the Laurent expansion around z = from both jo and ji , since the 
zq = contributions cancel in Eq. (|33|) anyway. 

Now, choose for F two straight hnesQ at angles it?, < < n/A, running 
to (from) the origin from (to) infinity, and rescale z . This gives 

IniT) = A / e-^^"/'h{uR)du, (40) 
Zm Jr 

with the same contour F in the u plane. 

The whole R dependence of the integrand is now in h{uR). When R oo, 
poles on the real axis condense into a cut. However above and below the 
cut the integrand simplifies considerably for large R. We exploit this by 
deforming the contour into two lines parallel to the real axis at distance e 
and a vertical section along the imaginary axis. Due to the Schwarz reflection 
principle, contribution from the vertical section vanish. Hence Ir{T) = + 

with 

4iT) = tR^^ f e-^(^±-)'/^Mi?(p ± te))dp, (41) 

At large R and fixed e trigonometric functions are dominated by the terms 
growing exponentially with R and in fact cancel in the ratios in Eq. (^91). 
Consequently we have simply (e.g., above the cut) 

z — 2i 

lim h(R(p + ie)) = — -, z = R(p + ie). (42) 

i?^oo, £>0 fixed ^ ^ " z{z — l) \-r / v / 

and symmetrically below the cut. 

Now we can move toward the upper and lower ends of the cut by taking 
e — 0. The imaginary contributions from upper and lower contours cancel 
and our main result reads 

/(r)= lim^ fe--"/'^-!^. (43) 

This can be easily calculated with the aid of the proper time integral repre- 
sentation 

I{T) = lim - / e-^"^"'^ / e-^(^'p'+^)rfsrfp (44) 
J. R e~^ds 1 ds 1 



^Contribution from the section of the big circle is negligible at large R. 
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which finally proves our conjecture based on Fig. [T| The limit singles out 
only the contribution from the p = state. In fact Eq. p3|) is nothing but 



im = / 6(p)e-^P'/^dp. (45) 
Jo 

This form clearly proves that supersymmetry is restored in the R —>■ oo limit. 
Densities of the bosonic and fermionic states are exactly equal for any nonzero 
energy, while the lowest state does not have a supersymmetric partner, which 
results in the point-like contribution from the zero momentum. This also 
confirms our earlier assertion that, in the D = 2 system, restricted Witten 
index is time independent even for the continuum spectrum [ED . 



Summarizing, supersymmetric quantum mechanical systems can now be 
solved with better precision in any dimensions. This provides a useful nu- 
merical tool for the nonperturbative solutions of matrix theories. Here it was 
applied to the N = 2 case. However, since the present recursive approach 
uses only gauge-invariant operators, it is even more useful for higher gauge 
groups. Finally, analytical calculation of the restricted Witten index for the 
D = 2 showed explicitly how the supersymmetry is restored while removing 
the infrared cutoff. 
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